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The BKZ elastie fluid theory is used to correlate experimental results obtained in biaxial strain 
and steady simple shear. With a heuristic potential function involving three material properties, 
excellent agreement is obtained between theory and experiment. In the special case where one of 
the material properties is dominant, the behavior in steady simple shear is calculated from dynamic 
measurements in the infinitesimal range and is compared with actual data. 
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1 . Introduction 

In a recent paper [l], 1 excellent agreement was 
shown between experimental results and predictions 
of the BKZ elastic fluid theory [2]. This theory in- 
volves a potential function, U, but leaves it unspeci- 
fied. For a given material, a knowledge of the results 
of a sufficient number of biaxial stress-relaxation ex- 
periments will enable one to predict with the BKZ 
theory the stress response to any other deformation 
history. However, if one knew a specific functional 
form by which U could be closely represented, then 
one would be able to correlate the behavior of differ- 
ent materials and different strain histories from the 
results of only a few experiments. 

Encouraged by these results [1], I constructed a 
form of V with which the BKZ elastic fluid theory 
could quantitatively describe biaxial strain at large 
deformation, biaxial creep, and simple extension of 
vulcanized rubbers. This form of U involves three 
material properties. In simple shear it can quanti- 
tatively predict non-Newtonian behavior, including 
normal stresses. The ratio of the shearing stress to 
the rate of shear depends on the rate of shear in such 
a way as to describe either shear thinning or shear 
thickening behavior or both, depending on the relative 
magnitude of the material properties. In the special 
case where one material property is dominant, one 
may use dynamic data taken at infinitesimal strains 
to predict the dependence of viscosity and normal 
stresses on rate of shear. This is presented in section 
4 of this paper and the agreement is excellent. 

I want to emphasize that the form of U presented 
here is heuristic. The purpose of this paper is to show 
that with a relatively simple form of U, one may use 
the BKZ elastic fluid to describe very well the behavior 
of materials which can be considered isotropic and 
incompressible. 



2. Theoretical Considerations 

The BKZ elastic fluid is a fluid with an elastic poten- 
tial. The effect of the configuration at time r < t on 
the stress at time t is equivalent to the effect of a stored 
elastic energy with the configuration at time r as the 
preferred configuration. The effect depends on the 
amount of time elapsed between r and t. The stress 
at time t is the sum of contributions from all past times. 
For an extensive description of the theory, we refer the 
reader to the initial papers [2] and [3]. 

A particular motion of the material may be specified 
in terms of the Cartesian coordinates xi of each particle 
at each time. Let X\,X2, Xs be the position coordinates 
of the particles in a reference configuration. Then, 
a motion is given by a set of functions 

Xi=Xi(X k9 1) U * = 1, 2, 3. (2.1) 

At time r, (2.1) becomes 

xi=Xi(Xk, t) i, fc=l, 2, 3. (2.2) 

If we eliminate X u X 2 , X 3 between (2.1) and (2.2) 
we may write 

Xi(t)=Xi(x k (T), t,r) i, A=l, 2, 3 

where xt(t) and %i(r) are the position coordinates at 
time t and r respectively of the same particle. The 
relative deformation gradients Xik(t, r) are defined by 



Xfk(t, t) i 



, Bxjt) 

dx k {r) 



Figures in brackets indicate the literature references at the end of this paper. 



The left Cauchy-Green tensor dj(t, r) is then 

Qj{t, T)=X lk {t, T)Xj k (t, t) 
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where repeated indices indicate summation over the From (2.3), (2.4), and (2.5) we get 
values 1, 2, 3. 

The principal invariants of Cij(t, r) are en \{t) — cr^(t) 



hit, t) - m& 

where £?=Xf(r, t) are the principal values of aj{t, t) 
and Xj is the stretch ratio in the xi direction. 

Assuming incompressibility, we have h(t, t)=1, 
and the constitutive equations for the BKZ elastic 
fluid become [4] 



xf(o yi(T)\i(r)-\idU , \%t) au- 
dit 



J^Lm(t) xf(f)xi( f )J La/, \1(t 



dr (2.6) 



and 

(Tr>{t) — cr 33 (r) 



< 



Xl(r) X?(t)X|(t)" 



A|(t) Xf(0Xl(f) 



af/ x|U)at/ 
La/. + X;(t) a/ 2 



dr. (2.7) 



0-y(f)=-pO(;+2/ 



a*/ . , at/ 



-. In the case of a single step stress-relaxation experi- 

~ l (t, t)\ dr (2.3) ment where X;(t) = 1 for times r smaller than zero and 
J ki(T) = \i(t) = \i for times t, r greater than zero, (2.6) 

yields 



where cr/j are the components of the stress tensor, 
p is a hydrostatic pressure, U is a function of A, h, 
and £ — t 

U=U(I 1 (t 9 T\I 2 (t 9 T), t-T), 

and c^ 1 ^, r) are the components of the inverse of the where 
matrix ||c y -(j, r)||. We may describe an isochoric ho- 
mogeneous biaxial strain history by writing for (2.1) 



<T\\(t) — CT; 



M- 



PME^S) m 



x?xi 



/; 



x,(f)=X 1 (r)^ 1 

* 2 (r)=x 2 ( c )A: 2 



* 3 W=X 3 ( f )* 3 = — Z 3 . 
AlA.2 

The matrix of the left Cauchy-Green tensor c//t, t) 
becomes 



and 



W=W(luh,t)=\ U(I u h,$(%, 



1 2 X?AJ 



(2.9) 



/ 2 = 



= — + — + X?X| 

x? xi ' - 



lko<r, r)|| = 



XjU) 
MM 



X|U) 
X?(r) 






Xf(T)XKr) 



XKOXKO 



Similarly (2.5) becomes 
CTjgQ)— cr 33 (t) 



^M^'f) « 



and 



/l(t,T) = 

/ 2 (f,T) = 



Xj(t) Xf(t) X^(t)XKt) 

X?(t) \|(t) Xf(t)XI( f ) 

Xj(7j Xf(T) X^(t)XKf) 

X*(r) k%t) X2(r)Xi(r)' 



^"XfXi 



Superficially, expressions (2.8) and (2.10) appear to 
be the same as the relations given by Rivlin and Saun- 
ders [5]. However W, here, depends on time as well 
as strain and is designed so that (2.10) gives the stress 
during stress relaxation. On the other hand, the W 

(2.4) of Rivlin and Saunders depends only on strain. 

In the case of simple shear, we introduce an orthog- 
onal set of coordinates as shown in figure 1, where %\ 

(2.5) * s tne direction of motion of a particle and x 2 is the 
direction of shear. In the case of a steady shearing 
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FIGURE 1. General coordinate axes for description of flow. 
motion with constant rate of shear, y, (2.1) reads 

Xi(t)=Xi(T) + (t-T)yX2(T) 

X2(t) = X>>{T) 
X:M = X:i{t) 
Thus, 



I Cij(t, t) I 



and 



\cTj l (t,r) 



l+y 2 (^-r) 2 y(t-r) 

y(t-T) 1 

1 

1 -y{t-r) 

-yO-r) l+y 2 (*-r) 2 

1 



vhere 



It(t,t-{)=h(t,t-& = 3 + y*e 



U = U(I,,h,{) 



. ., C»(dU dU] ., 

These relations hold independently of the form of U. 

3. Experimental Procedure 

From the theoretical considerations of the previous 
section, we see that we may determine dWjdU and 
dW/dh as functions of /i, 7 2 , and t from data taken in 
single step stress-relaxation experiments in biaxial 
strain. For vulcanized rubbers we may regard the 
long time isochrones in creep to within a good approxi- 
mation as isochrones of single step stress relaxation. 
This is true for a material only if the deformation at 
constant load remains almost constant for long times, 
although not necessarily to infinite time. For this 
reason we elected to do our experiments on vulcanized 
butyl rubber. The experiments were carried out 
using a test piece in the form of a square sheet having 
sides of 8 cm and a thickness of about 0.07 cm. The 
test piece was cut and marked in a fashion described 
by Rivlin and Saunders [5|. One square surface of 
the sheet was marked in ink with two sets of four 
parallel straight lines so as to form a square grid with 
1 cm spacings. In drawing the outer lines of the grid, 
great care was taken that they be straight and form a 
perfect square. In figure 2 is shown part of the ap- 



Substituting in equation (2.3) with t — T = f, one gets 




FIGURE 2. Schematic diagram for the biaxial extension experiments. 
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paratus with the sample S. As can be seen from the 
figure, the sample was stretched in the A— B and C— D 
directions with the aid of strings. One end of each 
string was tied to a lug of the sample and the other 
end was attached to a weight or to a turnbuckle. The 
strings were made from unspun fibers of Mylar. For 
the study of the variation of dW/dh and dWjdLi with 
respect to time, the two sets of the three middle lugs 
in the A and C directions were tied to strings support- 
ing weights. The other 14 lugs were tied to strings 
in which tension was controlled with the turnbuckles, 
which could be adjusted to keep the ink lines on the 
sample straight and parallel. With the aid of a two- 
way traveling microscope, whose axes of travel were 
set parallel to the stretch directions, we could check 
the uniformity of stretch and measure the extensions. 

In our other sets of biaxial experiments, we loaded 
the sample by first stretching it a predetermined 
amount in the D— C direction and then applying weights 
in the A-B direction. This was done in order to avoid 
difficulties due to the history of loading. In these 
experiments, only the three middle lugs in the A direc- 
tion carried supporting weights. The strings attached 
to the other lugs were adjusted to keep the lines 
straight and parallel to the stretch directions. The 
final readings were taken after 18 hr from the time of 
loading. After each measurement the material was 
allowed to relax for 24 to 48 hr before another loading 
was started. Thus, the values that we obtained can 
be considered as isochronal values of stress relaxation 
at 18 hr. 

In order to be able to compare with experiment the 
predictions of our theory for simple shear flows, we 
performed dynamic and constant rate of shear meas- 
urements on solutions of polyisobutylene B-140 in 
Mentor 28 oil. Two concentrations were used. They 
shall be designated as 10 percent and 5 percent. We 
do not know the actual concentration accurately, be- 
cause we lost an unknown amount of solvent while 
preparing what was to be the 10 percent solution. 
We do know that the ratio of the two concentrations 
is two to one. The dynamic data at very small defor- 
mations were obtained through the cooperation of 
R. W. Penn, using a torsion pendulum at the W. R. 
Grace Laboratories. The torsion pendulum is essen- 
tially the same as the one described by Morrison, 
Zapas, and DeWitt [6]. The data on viscosity as a 
function of rate of shear were obtained in a capillary 
viscometer. 



4. Experimental Results and Discussion 

The purpose of this paper is to show that with a 
relatively simple form of the potential function U, 
the BKZ elastic fluid can be used to correlate different 
types of behavior of elastomeric materials. The 
heuristic form of the potential function U which I 
shall use here involves three material properties ot(t) 9 
(3(t), and c(t). These material properties are positive 
monotonically decreasing functions of time. The 
form of U is given by the following expression: 



~U=y (/,-3) 2 +4.5/3' In ( /|+ ^ 2 + 3 



(iTTEJ+c'Vi- 3 ) ( 4J ) 



where 



, da(t) d(3(t) , dc(t) 

a = — — , b = — — and c =— — 
dt dt dt 



From equations (2.9) and (4.1) we get 

dW - n , Y , 4.5)8 24Q3-c) , 

dW = 4.5/3 24(/3-c) 
dh /i+/ 2 + 3 / 2 +15 



(4.2) 



(4.3) 



where it is understood W \ a, /3, and c are functions of 
time. 

In a pure shear single step stress-relaxation experi- 
ment with \ 2 = 1, and /i = / 2 , we get from eq (2.8) 

0-11(0-0-33(0 ^ IbW dW\ 
X2 _l \Bh dh) 

1 M 



2 a(/i-3) + 



9/3 



+ c (4.4) 



2A + 3 



Since a(t) is taken to be a positive monotonically de- 
creasing function of time, or zero, eq (4.4) says that 
•r 1 . 0-11(0-0-33 (0 1 

it one plots at constant t versus 



one should get either a curve which is concave up- 
ward or a straight line. In the case of a straight line 
the slope is equal to 18/3, the intercept, 2c, and a 
equals zero. 

In figure 3 we show the data of Rivlin and Saunders 
[5] on pure shear for vulcanized natural rubber. In 

this figure we plotted —y ~ — a(I t — 3) versus 



2 Af- 



K 



2/x + 3 



for a = and a = 0.06. As can be seen for the case 
where a = 0, the curve is concave upwards. Here 
a = 0.06 was found by trial and error. Actually, if 
there is a well defined minimum, one could obtain a 
by following the procedure presented in a previous 
paper [3] . So it is evident that by using the relation 
(4.4) one could get the three material properties from 
pure shear experiments. 

In figure 4 we present a similar check of the ade- 
quacy of the assumed form of U by plotting 



0-11(0 — 0-33(0 



"'~w) 



-a{I\ — 3) versus 



1 + X| 

/1+/2 + 3 
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FIGURE 3. Pure shear data ofRivlin and Saunders on vulcanized 
natural rubber. 

Open circles, a = 0. The abscissa is given in kg per cm 2 . 



Figure 4. Biaxial extension data of vulcanized butyl rubber. 

The abscissa is given in (dynes/cm 2 )10~ 6 . 



using data obtained on vulcanized butyl rubber in 
biaxial and simple extension deformations as de- 
scribed in the previous section. For butyl rubber, the 
value of a(t) is small and estimated to be 0.015. Since 
fi and c are equal, as determined from the pure shear 
data, in a plot of this type the experimental points 
should fall in a straight line. Considering experimen- 
tal difficulties and uncertainties, the agreement is 
excellent. 



TABLE 1. Biaxial creep of vulcanized butyl rubber 



Time 


/, 


I> 


f*- 


ff~ X 10 ~ 6 


Hours 






Dynes/cm 2 


Dynes/cm 2 


3 


3.1043 


3.1050 


2.15 


-0.10 


20 


3.1125 


3.1125 


1.91 


.07 


23 


3.1134 


3.1131 


1.83 


.14 


47 


3.1160 


3.1156 


1.80 


.16 


66 


3.1177 


3.1172 


1.78 


.17 


90 


3.1199 


3.1191 


1.74 


.19 


117 


3.1211 


3.1203 


1.71 


.20 


164 


3.1233 


3.1224 


1.69 


.21 



In the course of our experiments in biaxial deforma- 
tion, in a different set of measurements we observed 
a negative value of dW/dh at very small extensions, 
while at higher extensions dWjdh was positive. This 
was observed in three samples with different degrees 
of vulcanization. For a further study of this peculi- 
arity we selected a relatively high cross-linked speci- 
men of butyl rubber and we studied its biaxial creep 
behavior at small deformations. In table 1 we show 
the calculated values of dW/dh and dW\dI 2 as func- 
tions of /i, 7 2 , and time. This table shows large varia- 
tions in dWldI 2 for small changes in h and / 2 . A plot 
of dWIdh + dWjdlz versus the logarithm of time is 
shown in figure 5. We see that even in greatly ex- 
panded scale for dWIdh+dWjdh the points fall in a 
straight line. Moreover, after 150 hr the data still 



Log [t, hours! 



FIGURE 5. Biaxial creep data on vulcanized butyl rubber. 



did not indicate any suggestions of leveling off. This 
shows that we cannot neglect the behavior with re- 
spect to time. From single step stress-relaxation 
experiments in simple extension in vulvanized rubbers, 
we observed that (3 decays with time much faster than 
c. From (4.3) we can see that dW/dI 2 will be negative 
when /3 is larger than 1.6c. One can see, at least 
qualitatively, that what seemed to be a paradoxical 
inversion in sign for dW/dI 2 is predictable from eq 
(4.3). 



4.1. Steady Simple Shearing Flow 

We can substitute eq (4.1) into eq (2.11) to get 

o- 12 =- 2 ( X laWe + %^§- + c'(fll y?d£. 

Jo I 9 + 2yt J (45) 
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Since the viscosity rj(y) is defined to be the ratio of 4.2. Comparison of Steady Simple Shearing Flow 
<Ti2 to y, eq (4.5) can be written 

i?(y) = — 2 



a '^r^ 2 + 9^^+^^) 



With Measurements at Infinitesimal Deformations 

For a simple shear deformation, which is specified 
in terms of a single parameter, say y=y(r), (2.2) may 
be written 



(4.6) 



(4.8) 



By considering special cases for which the dominant 
term is that containing a'(£), /3'(f), or c'(£), one can 
see that (4.6) could predict a viscosity independent of 
rate of shear, or viscosities as observed in shear thin- 
ning or shear thickening materials. 

In figure 6 we show schematically the type of be- „ 
havior predicted from eq (4.6). For the sake of To represent single step stress relaxation,^ we take 



x l (r)=X l ^y(r)X 2 

#2(t)=X> 

x 3 (t)=X 3 . 



simplicity we can write 



(4.7) 



y(r) = for r < and y(r) = y = constant for r > 0, 
and we get 



(4.9) 



r,dy)=-2 a'(£yy*e<% 



w) 



y)=-2 f" 

Jo 



1+| y^ 2 



« 



Vc=-2 \ c' 
Jo 



£*<# 



(4.7a) 
(4.7b) 

(4.7c) 



1 + y 2 


7 





y 


1 











1 


l 


-y 





-y 


1 + y 2 





o 





1 



where I\ = I>> = 3 + y 2 
and (4.9), we obtain 



o-i2(0 = 2y 



From (2.3), (2.9), (4.2), (4.3), 



*"W + ^ + *) 



(4.10) 



We observe that in the case where T7a(y) — T?/3(y) — 

one gets a viscosity independent of rate of shear. The ^ of ^f) for vanishi ives us the rdaxa . 

When only 7] a (y) = 0, one gets a behavior shown in 7 

curve II of figure 6. Curves III and IV represent the tion function, G(t), for infinitesimal deformations: 

cases where T)«(y), ^My), and r) c all contribute to the 



viscosity. Naturally, curve IV shows the case where 
7] a (y) is the dominant quantity. 



G(t) = 2p(t) + 2c(t). 



(4.11) 




From the general relation of linear viscoelasticity [7] 
between periodic and steady-state functions 

G'(o>) + ia)r)'(co) = ico J Ofye-^dt, 

we can express r)'{(o) in terms of fi{t) and c(t) as 

17 '(w) = {2/5(0 + 2c(£)} cos cwfcfc. (4.12) 

Jo 

Integrating (4.6) by parts we get 
W=2 f 3C fo-^--.^. 9j8(f)(9-2yy 



T?(y) = 



3a(£)y 2 £ 2 + : 



(9 + 2y 2 £ 2 



(4.13) 



LOG 



y 



Several interesting observations can be drawn re- 
garding the behavior of a material which can be de- 

FlGURE 6. Schematic representation of forms of steady shear Scribed by a potential function U of the form of (4.1) 
viscosity curves which can be predicted by eq {4.6). by comparing (4.12) and (4.13). First, it IS obvious 
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that in general 77(7) cannot be predicted from measure- 
ments of T7'(co), since the term in a(t) does not show 
up in r)'{(o) at all. Also, /3(t) and c(t) cannot be sepa- 
rated by measurements of r) f (a>). However, if we 
encounter a material for which either /3 or c is the 
dominant term in (4.13) and (4.9) then in principle 
77(7) can be calculated from 77 '((o). If c is the dominant 
term, we will have a situation in which 77(7) is inde- 
pendent of rate of shear, but 77 ' (co) may vary with fre- 
quency and normal stress effects may be observed in 
steady shearing flow. If /3 is dominant, both iq'(o)) 
and 77(7) will vary with a> or 7. 

For many materials a(t) is negligible, but ordinarily 
both P(t) and c{t) contribute to 77(7), with fi(t) contribut- 
ing the dominant term within experimentally acces- 
sible rates of shear. However, at very high rates of 
shear, any nonzero c must become dominant. For 
the range of rates of shear for which 



/; 



9^)(9-2y^ f» 

(9 + 2y 2 a 2 * Jo *'* 



we may utilize 77 '(co) to evaluate /3(t) and calculate a 
lower bound for the measured 77(7). The calculated 
77(7) should be in close agreement with that measured 
at low values of 7, but would fall below the measured 
values at high rates of shear. 

For the actual comparison of the two measured func- 
tions 77' (co) and 77(7), it is better to formulate our ex- 
pressions in terms of relaxation spectra corresponding 
to fi(t) and c(t) entirely analogous to the spectrum repre- 
sentation of G(t). The relaxation spectrum F(t) may 
be defined by [7]: 



G(t)= f ' 

Jo 



Fi^e-^dr 



(4.14) 



and may be expressed as the sum of two terms F^r) 
and F c {t) defined by 

2/3(0=[ ^ F^e-^dr 

2c(0= f * F c (r)e-*dT. 

Jo 

In terms of (4.14) we have 

r -.{^(r) + F c (r)}w 2 T 2 rfr 



CM 



Jo 



-hcuV 
r)}rdr 



(4.15) 



vM=r {F ^:^? }7rfr (4.16) 

Jo H-co-T- 



and in the case where a(t) = from eq (4.6) 
f f°° xe~ x dx 



v(y)~- 



**j. i + |*w 



+ tF c (t) dr. 




3 



Figure 7. 
28 oil. 



Table 2. 



Dynamic data on 5 percent solution o/B-140 in Mentor 

Solid circles represent the values calculated from data given in table 2. 



Relaxation spectrum for 5 percent solution of polyiso- 
butylene B-140 in Mentor 28 at 25 °C 



T 


H(r) 


19.0 


0.0 


15.2 


.1 


10 


.5 


3 


7.3 


1.0 


39 


0.3 


120 


.1 


245 


.03 


410 


.01 


420 



We calculated H(t) = tF(t) from measurements of 
T7'(o)) and G\o)) on the five percent solution of B-140 
by an iterative method which will be described in 
another paper. In figure 7 we show the dynamic 
rigidity and viscosity as a function of frequency. 
The black points represent points calculated from 
H{r) as obtained by our iterative method and tabu- 
lated in table 2. The agreement indicates that we 
have a good representation within the range of meas- 
urements. In figure 8 we show the viscosity at steady 
shearing flow as a function of rate of shear, with the 
open circles representing values calculated from eq 
(4.17) assuming that the contribution of the integral 
of F c (r)r is negligible. The agreement between ex- 
perimental and calculated values is excellent. 

The same arguments can be used for the deter- 
mination of the normal stress differences cru — a, 
o"22 — cr, and a 3 3 — cr (where 3o-=<x n + a 22+0-33) as can 
be seen from eqs (2.12) to (2.14). An interesting result 



is that the limiting value of the ratio 



0*2 



-0-33 



at small 



(4.17) 



CT\ 1 — O-33 

rates of shear is 0.46 in the case where /3 is the domi- 
nant term. This compares very well with the value 
of 0.4 reported by Markovitz [8] for a 5.39 percent 
solution of polyisobutylene in cetane. 



531 




v(y) 



-2 



I 

LOG |y=Gu,sec H ] 



Figure 8. Dynamic viscosity and steady shear viscosity data on 
5 percent solution of B- 140 in Mentor 28 oil. 

Solid circles are the calculated values. 

In conclusion, we should emphasize that all these 
derivations were obtained with the assumption of 
incompressibility. In reality dW/dli and dW/dh de- 
pend also on 7 3 , the influence of which can not be 
evaluated from the experiments reported above. 
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